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I—ON THE RESIDUAL CALCULUS. 


Caucny, in several articles in his “ F/xercices,” has developed a 
new Calculus, to which he has given the name of the ‘ Calcul des 
Résidus.” Few persons appear to have followed the author in 
the use of it; and in our own language I do not recollect having 
met even with an allusion to it: yet it deserves some attention at 
the hands of mathematicians, were it only from the deserved ce- 
lebrity of its author. I propose to give here merely a slight sketch 
of its principles, with a few examples ; and if any one should be in- 
duced from this specimen to wish for a more complete knowledge 
of this Calculus, I must refer him to the Memoirs of its author. 

The Residual Calculus bears a certain analogy to the Differential 
Calculus, for in both the object is to investigate the nature and 
properties of certain functions which appear in an indeterminate 
form, but yet have finite values. 


Let f(x) be a function of 2, which becomes infinite when # =a; 
that is, let a be a root of the equation 
1 
(ys 6; or a> = 0 Cl). 
I) 7@) 

Then the expression (a —a) f(x) becomes indeterminate when 
«=a, as it takes the form 0 x @, though it may still have really 
a finite value. Suppose this to be the case, and that it is a certain 
function of a which we may call R(a@). Then we have 

(z — a) f(x) = R(x) when zx =@............ (2), 
R (a) 


which gives f(z)= ert WHO DoD siiccicicess LO) 


U 
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This function R (a) is called the residue of f(x) with respect to 
a. The operation of finding the residue is called the extraction of 
the residue, and is represented by the characteristic symbol ¢, so 
that we write 

R(a) =¢ f(2). 
In equation (3) let us suppose x =a+A. Then 


R(a + h) 


] 
flat hy = i = i sR (a) +h R'(a) + &e.3, 
t 


when xz =a, or A= 0. 
Whence we find the residue R(a) of f(a) is the coefficient of j 
: h 


in the expansion of f(a-+h), where a is a root of the equation 
S(zj=o@. 

We have spoken only of one root of this equation, but there 
may be several, and for each of these there is a corresponding 
residue. The sum of the residues corresponding to each root of 
the equation is called the integral residue of f(x). To distinguish 
the integral from the partial residues, I shall suffix the root to the 
Residual Symbol in the latter case.* Thus, ¢ f(«) will represent 
the integral residue of f(x), g. f(x), ¢,f(2), ¢.f(@), &e. will 
represent the partial residues with respect to a, b, ec, &e. these 
being roots of the equation f(r) =o. 


Let us suppose, now, that instead of the expression 
(x — a) f(x) 
having a finite value when x = a, that it is the expression 
j m a 
(x — a)"f(x) 
which has a finite value when =a; in which case f(v) = @ is 


said to have m equal roots. If, as before, we suppose R(a) to be 
the finite value, we shall have 


f(x) = R (2) when r= a, 
(2 —a)™ 
or putting 7 = a-+h, and expanding R(a +h) by Taylor's 


theorem, 


hp Ai OR™D (a) 
— i xe, — = 
f(a+h) a{R(a) + i (a) + &e. + 1.2...(m — 1) + se} 


. ee 

The coefficient of , this expansion is 
R™~!)(a) 

12. .. (m a i) ‘ 


which is the residue of the function f(a) corresponding to 2 = a. 





* This is not Cauchy’s notation. 
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Since R(a) = — = hm fia + h), 


R=) (a) = sh" f(a rN hy when h = 0; 


“ipa I 
and therefore we have 
: : 1 — 1 
(4) gf (x) = — (m— 1) an He sh” f(a+ h) when A= 0. 


I shall now ae to show how the extraction of residues may 
be simplified when f(«) has particular forms: but previously it is 
necessary to explain the following notation. When f(2) consists of 
distinct factors, which are functions of a, the roots of f(x) = @ 
may be the roots of any one factor equated to «. 

When then we wish to indicate that the extraction of residues 
has reference only to the roots of one of the factors, that factor is 
enclosed in brackets [ ]. As, for instance, if f(a) = @ (x). (2), 
and we wish to represent the extraction of the residues with re- 
ference to the roots of g(a) = ©, we shall write 


§ [o(x)].Y (2). 


] 
If a factor is of the form oa)’ we should properly represent 
av 


the extraction of the residues with respect to it by 


r es 
Lawl , 


But since 7) = * gives the same roots as @ (a) = 0, which, 


therefore, is the equation to which we look, we shall find it con- 
venient to use the notation 
¥ (2) 


[y(x)]} 

Hence, if a be a root of f(a) = #, the residue with regard to it 

(x —a)f(a ’) 
|v — a} 


€ eae) =hAf(a+h) when h=0, 


if there be only one root of f(a) = 2; and 


( r—a). f(: i= 1 qn 
[x—a ] 1.2...(m—1) dh" 


if there be m roots equal to a. 


¢ 


may be expressed by § instead of ¢, f(a), so that 


sh f(a +h) when h=0, 





N (x) to be of the fi (2) 1 that a i t 
rs se f(x) to be 1c form ——~ , and that a is a r 
Now suppose f(a) to be o rm F(z) and that ¢ 4 TOO 
of F (x) = 0, then, as F (a) = 0, we have 
F(ath)=hF(a4 th), 
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and therefore 
At(a+h) 


“S(4 +) = TE G+ Oh)’ 


So that when 4 = 0, 


ro f(x) a f(a) 5) 
- [F(2)] F’ (a) tesseeeesees CD) 
As an example, take souks whence 
r—a 
am . ” 
cs. fe—a a a} = a”. 
i ,@—pr+gq ” a> — pa + q 
; ** (a? — a®] Ya : 
; am as a” 
a [(a — a)(x — b)) ” ian 
f(a) 


If the fraction be of the form » wher: 


(« — a)(#—b)...(@—71) 
there are r factors in the denominator, we may extract the residue 
with respect to each factor, and the sum of all these partial residues 
will be the integral residue with respect to the denominator. 
Hence 


f(x) : f(a) 
9 {(a—a) (a — b)...(x — r)] (a — b)(a—c)...(a — 71) 
f(b : f(r 
+ (b it + &. + (r— a) (7 ‘ = ro q) ; 
When the denominator has m roots, each equal to a, 
f(v) an 
‘5 [F(ax)} ~ 4.2...(m —1) dh™ 
when 2 = 0. 
f(a +h) 
F(a + 4) 
a © Lai Oa 
1.2...(m—1) 1.2...m ; 
and as these are m equal roots of F(a) = 0, we have 
F (a) =0, F (a) =0.... F"™—)(a)=0. 


c.Siz) = shm f(a 4 hype 


But f(a + h) 


» and 


| 
F(a+h)=F(a) + F(a); +8. + 


hm 
Hence F(a +h) = : SF™ (a) + Fi) (a) 


= . &e. 
1.2...98 $? 


h 
m+) . 
7 m on dm 3 
or F(a + h) - fii qa* F(a 4 h);* 


ne 


/ 
f(a +h), 1 mean to express the value which 
da™* 


* By the notation 
qm 


ia f(x +h) takes when a a 
TX 
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whence we find 
ry f(x) Pe ] d=! 
[F (x) ] 1.2...(m— ] ) dhn i, ‘2 


hn f(a+th) 








“dh = F (a+h) 
when A =0 ....... cc cee ace (6). 
When I (a+) is of the form (@—a)", this expression is con- 
siderably simplified. For 
d™ 


= (2 ~ a)" = m(m — 1)... 2.1, 


m dm : 
and [diem F(a +h) = hk”, so that 


da” 


on ore ee 7 
“4 (¢—a)™— 1.2...(m—1) dh™-! ¢ (4 + 4)3 when A=0...... (7). 
ax” n(n — 1) (2 -— 2) 


~ 3 
(# — a)! 1.2.3 , 


a* 


_ 
rhus §, 


I shall now proceed to the proof of one of the fundamental 
theorems of the Residual Calculus. It is, that 
f(z) _ 


S (x) -€ a re (8), 


Y(«) being a function of 2, which takes generally a finite value 
when 
t=, = My» a, &e., 
(,, A, a,, being the roots of the equation f(a) =~. 
We have generally 
y R (a 
f(z)= : when «= 4,. 
(x oa a,) 
Expanding R(x) = R(a, +A) by Taylor's theorem, and putting 
x—da, for h, 


1 
R (a,) l R’ (a) 


; r) = &C. 
t( ) (x a a,)” l (x — a)” 1 + : 
] R(™—) (a,) 
? | d 
t 1.2..(m—1) a—a, F EE) ore aeeee (9)s 

where y (a) is a function of 2 which generally takes for a = a, the 
fil R(™ (a,) 

finite value Ta 


But we have also from the definition of a residue 


R™ (a,) Pe R (2) 


(10), 


1.2 ...90 (e—a,ye*) 
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Hence 


and similarly for the others. 


l R(z R (2 
S(2) = is =a)" 8a, : ¢ 2 + £., (2 “ ¥ (a — a) 
; R (z) . : 
+ $a (2—a pe —ayP+4 &e} 
R(z) 


=a (7 -a,)m za, ni? ~a,)™"! 4. (2—a,)™-*(a—a,) + &e.3 +(x) 


R (z) {(@ —a,)™—(z—a,)" oe 
=$u \(a—a, "(z—a, ym | . = a: free) 
R (z) . R(z) 


=fa(z. —z)(z—a,)™ (x—a,)" “zz t v(z). 


Now the second term of this expression is equal to 0, because 


R (z) 
= does not become infinite for z=a, and consequently has 
no corresponding residue. We have therefore 


Re ‘(z 
since R(z) =(z — ~ am fiz). 

Hence, in order to deduce from f(x), which becomes infinite for 
x =a, a function of 2, which shall remain finite under the same 
circumstances, we must subtract from f (x) the partial residue with 





SC Te ee . " ee . 
respect to a, of f ) . Similarly, if we wish to find a function of x 


which shall a hesuns infinite when a =a,, =a,, &e we must 


these roots. So that if we make 


"pede a 


f(x) — = 9(2), 
@(«) will have a finite value when # =a,, =a,, =a,, &c. 
Zi (z)] . 
The notation § L. oF it must be remembered, means the sum 


of all the residues with respect to the roots of the equation f(z) = « 
From this theorem we can deduce a number of important re sults. 
If (2) is 0 fraction, as 
Jf (2) is a fraction, a: F@) 
kind, the denominator of which must never become 0, since the 
fraction ¢ (a) can never become infinite. 


, 9(2) must be a fraction of the same 


Therefore the denominator must be independent of 2 or constant, 
and @(a) an integral function of 2. If F (v) be of higher dimen- 
sions than f(a), f(x) will vanish for 7 = ©, as likewise will the 
residue, so that we must also have o(a2)= 0. Hence in this case 

LA (z)] 
f(r)=§ . i if), 








subtract from f(a) the partial residues with respect to each of 
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On the Residual Caleulus. 


which gives the means of decomposing a rational fraction unto sim- 
ple fractions. 


x? — Tax + 1) 








As an example, take the fraction "= E=2)@—8) = f(r). 
. —7z+1 
f(x) =S7, —z) Teer —2) (z—3)]_ 
—7z ve ] 22— 72+] 
¢ + $2 , 
EEG ha al 
2 as 7z }- ] 
ts @ ae NE- HEH) 
5 ] «..* .. # 1 
8a) sie x—2 224-3 
Again, take ——— _—_- f(x) 
ss (w—1)?(a4+1) ° 
l l 
fle) = So (=a e—ipP etl) + & ez) (e—1P (4 
l d 1 
= aa. aes vhen h = O 
“t-77* @6-i-baaa 
l ] 1 


1 ee ee 
4e+1 +3 (e—1)? *a—1 


. re f (a) i i ial Maia 
When the fraction is of the form a» We can easily deduce 





(x—a) 
by this method the series of component simple fractions. 
— 
or (x sane a)" — @a (a Av z) (z nae a)" 
im—1 
= ie al ( =a h) ) when h = 0; 
1.2...(m —1) dhk™!'\xa—a—h 


and if we actually ptm the operations indicated, we get 





f (a) mee Boe — ~ fim?) (a) 
gaaP Tae fea) Yeap 
(m3) (q 
+ (m— 1) (m— ¢ ee &e. 
f " 
+ (m — 1) (m — 2) ...2.1 ana! ectaniedy (12). 


This also appears from (9) by inverting the order of the terms. 


By the application of the theorem (11) we can easily prove La- 
grange’s formula for interpolation. For if 
f (a) 


(2 — 2,) (w@ — a,) ... (w — 2,,) 





f(z) = 


’ 
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we have 
e f (z) ] 
* (z — 2,) (2 — a)... (2 t,) &2— 2 
. f(z) | 
os (2 z,) (f — a) ... (8 x.) 2#—2 
f(z) l 
s®o - ~— & 
+ ~(z — 2,) (z — 2,) ... (¢ Ln) L— 2 lies 
,; f(z) 1 
“™ (z — 2) (z — 2) ...(2 — 2) 2— 2 
f (2,) | 
oa - &e. 
(2, — 2,).... (2, Lt fy t 
a f(2,,) 1 
(¥,, — By) oo Bey Ln, © — Ly, 


(x — x,) ... (wv — 2,,) 
(=, ~ &) ... (a, — 2.) 


4 (x7 —a2,)..%— a 


whence f(x) = f(a) + &e. 


‘m—| f(a - 
(x, it «,) one Ly ak . 


- F : m=) 
Again, taking the equation 
t(e)= ¢ OQ), 
r—2z 


and multiplying both sides by 2, it becomes 


, Lf (2)] 
zf(z)= —— 
Secs 
x 
On making «= ©, f(x) = 0, as the denominator of f(2) is 
supposed to be of higher dimensions than the numerator ; therefore 
x f(x) may have a finite value. Let this be V, then 
SIZ) = V 20. cccsceees (13). 
If V = 0, which will be the case when the degree of the de- 
nominator of f(a) surpasses that of the numerator by more than 
unity, then we have simply 


GS(Z) HO ceceeeees oe (14) 
Let S(#) wi (a - x) (a _ 2) _ (x “ a,,)? 
we find that 





a," = —— ay — + &e 
(a, — 2,) (a, — t,,) (t% — 2,) ... (4, — 2,,) 
n 
x 
E — m == 1 or = @, 
(z,, aed x,) eee (z,, = Ton) 


according as nm =m— 1, or as m < m— 1; since in the former 














Un ) : 


Qs 0) a) 
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case we have x f(a) = 1 when a = «&, and in the latter « f(x) = 0 
when 7 = @. 


We have as yet limited ourselves to the consideration of proper 
fractions, but the Residual Calculus may be extended to all others. 
To do this, we must first prove the following Theorem, 


f) 


Ef(z Jab Tey os eee AB). 
We know by (8) that 
] ] 
BA f\= 
‘oe 





uC «d's*]} (uw — 8) 


where U is a function of «, which remains finite for « = 0. 
fk 
Let there be m roots of f (=) =o, each equal to 0, and let us 
oN 


] 
expand ———; then 
u—s 


+U. 


aig e+ 8 


“a a? Ts*] “Tey 7] + Set ome & [s*] 


Multiplying by wu’, substituting z for - , and making Uw? = ¢(z), 


Ha) fia) mss) 


&e am ig - 
ay + Spey tee + 8 
Now, extracting the residues with regard to z, and observing that 
¢ (+) =1 g1=0 gs=0...g2s*=0, 
1 
(3) 
THT 


Now ¢(z) = Uw?, and U has a finite value for «= 0, there- 
fore Uw = z@(z) must vanish. Hence by (14) 


¢ 9(z) =z ¢(z) = 0. 


Hence, finally, we have 


i) 





f(z) =28 +9(2). 


we find of (z) = § —— + £¢(2) «-- «. (16). 





ef(z) = 8 


[s*]  * [2*] 
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Si) 


In this last formula substitute = 





for f(z). As the residue 
[Se) 


=o, and one of 


is taken with regard to all the roots of —-— 


these is # =z, we may separate that suid residue from the 
others. Thus we have 


pf 26 JO 4g LOI 





*n 
~ 


nee 


i) 
(f(z So) Fries) GTO==) 

1 
or f(x) = g LAC) Ig a olla (18), 


which is the equation we ought to substitute for the formula (11), 
if f(z) becomes infinite when z = 0. 


=f(z) + § =—— 


so that, by the last dinsenies 


S(@) +85 





If f(a) is an integral function of x, the residue § (ey) vanishes, 


since f(z) has a finite value for any value of z, so that the expres- 


sion is reduced to 
1 
f(z) 


2) = nee 9). 
S(2) * — 2x) ( 

If f(x) be of the toe 5 ie t and F(a) be of higher dimensions 
than f(x), the expression (18) is, as was proved before, reduced to 


$ Ls 1, which enables us to decompose a rational fraction. 


But if F(z) be of lower dimensions than f(a), the function ray 
is, by (18), divided into two parts, of which § LA) is the sum 
of the rational fractions remaining from the division of f(a) by 


AG) 


F(z), and § (sTC—se) is the quotient arising from that division. 
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Let us take, as an example, f(a) = a) 
cs (1 +24) a 1+-24 
f(®) = [z(1+2%)] (a—z) +& [2] (142?) (1 —z2)’ 


Now 








| Ay ee | a 
Sx 1+22)](a—z) 61 $2) (waz) iis li +2*)(a—z) 
; 14-24 
+8-— 1 aF2)@=2) 
— a 
2 r+— hy ee 1 
1424 d (14h) 


and § when h=0, =2.. 





[2*] (1+2*)(l—zx) ~ dh (1+4*) (1—ah) 
Therefore, finally, 





1+2* 1 1 1 
9 =wv — — 
x(1+2*) +; r+—?] r-—'*] 
] 2x 
ail ae 


The length to which this article has extended, renders it necessary 
to postpone to a future Number other illustrations of the use of the 
Residual Calculus. 

D. F. G. 


IIl.—_TRANSFORMATION OF CERTAIN ANALYTICAL 
EXPRESSIONS. 


WE frequently require in Analysis to find the sum of the squares 
of three quantities of the form 


ay — bx, bz — cy, cx — az, 
and a slight artifice enables us to do this very readily. 
By adding and subtracting a*a* + b*y* + c*z*, the sum of the 
squares becomes 
(a® 4-6? 4c?) (a? + y* +. 2*)—(ax+ by--ez)?......(1) 
ax + by 4+ ¢z)? 
nee ee ae ( mE apy + 7} 
ax+-by +cz 
(@+PFEN (ype 


the angle between two lines which make angles with the coordinate 





is the cosine of 





Now the expression 
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axes, whose cosines are proportional to abc, «yz respectively. If 
this angle be called 6, the expression becomes 
a® + b* 4 c*) (2* 4+ y*® 4 2°) sin’ @......(3). 


We shall give some examples of the use of this formula. 








1. In transforming expressions from one set of rectangular co- 
ordinates to another, making with the original axes angles whose 


cosines are aa'a", bb'b’, ec'c", we obtain the relations 


bc! —c'b" ; ch" — be! - bc! — ch! 
i: miei edeeias wie 
and we wish to discover the value of D. Squaring these equations, 
adding and observing that a? + a’® +a’? = 1, we find 
8 § 
D? = (b'c! —c'b")? + (cb" — be")? + (be — eb. 
By formula (3) this becomes 
D? = (0? 40 4b") (2 4. c'2 + ¢!2) sin? 6, 
Now sin 6 = 1, since the lines are at right angles to each other, 
the new coordinates being rectangular; also b* + 6% + 6% = 1, 
and c? +c? 4 c’? = 1, so that D = 1, and 


U 
a = Ue’ — cb", a’ = cb" — be’, a! = be — cb. 






















2. Let a material point be acted on by any forces whose resolved 
parts parallel to rectangular axes are X, Y, Z; let p be the radius of 
curvature of its path at any point ayz, then we know that 
1 /ds\3 __ 

7 (at) = 
{ (5 d’y dy =) 7 (2 dz dzdy\2 /dzd*x dx d*z\2)3 
dt d2 dt dt de Ea) & ii a ZA) } 


= , dx , dy\* , dy . dz\* , dz , dx\*)4 

={(¥5 — xi) + (2% -~Y¥F) + (xz — 23) }. 
Therefore, by formula (3) 

1 /ds\8 a 

~(5) — VX24-Y?47? asin (ang. between tang. and resultant), 

or if R be the resultant, V the velocity, 

v* = p R cos(angle between resultant and rad. of curvature). 


Calling this angle 9, R cos ¢ is the part of the resultant resolved 


perpendicularly to the path, which is opposed and equal to the 
, er v 
centrifugal, which is therefore equal to—. p cos @ is the 4 chord 


of curvature lying in the direction of the resultant, and therefore 
v? = 2R 4 chord of curvature along resultant, 


or the velocity in any curve at any point is equal to that acquired by 
falling through } of the chord of curvature under the influence of 
the resultant force at the point, considered constant. 








If 


co- 
hose 
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3. If p, g, 7 be the angular velocities of a body round the coordi- 

nate axes, the velocities of a point (yz) parallel to the axes are 
gz — TY, TX — pz, py — qa. 
A point in the line whose equations are 

coy 2 

-- @.* 
has no velocity. This line therefore is the axis of instantaneous 
rotation, and p, g, 7 are respectively proportional to the cosines of 
the angles it makes with the coordinate axes. 

The velocity of any other point is the square root of the sum 
of the squares of the velocities parallel to the axes. If this be v 
we have by (3) 
v= Vp4g4r Var 4 y? 422. sin (angle between inst. axis and 

line through 2yz) 


And Vp?+@47%, being equal to the linear velocity divided 
by the distance of the point from the instantaneous axis, is the 
angular velocity about the axis. 

4. The expression for the radius of curvature of a curve of 
double curvature is 
s? $( Ady — Bdx)? + (Cdx— Adz)? + (Bdz—Cdy)?34 


~ A24+ B24 C? 





where 
A = dy @z—dz d*y, B = dz @x—dxd’z, C = dx d’y—dy dx, 
which are proportional to the cosines of the angles which a normal to 
the osculating plane makes with the coordinate axes. The nume- 
rator is, by (3), equal to 

ds? (A? + B? 4 C*)3 (da? + dy* + dz) sin? 0, 


when @ is the angle between a normal to the osculating plane and the 


tangent, and therefore is a right angle, so that sin @ = 1; and, 
as dx? + dy? + dz? = ds*, we have 
ds 


P(A? 4+ B24C2)4* 
A.S. 


IIl.—DEMONSTRATIONS OF SOME PROPERTIES 
OF A TRIANGLE. 


l'une following properties of a triangle are not perhaps generally 
known. If we join D, E, F, (Fig 1,) the points in which perpen- 
diculars on the sides from the opposite angles intersect the sides, 
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the triangles DEC, DBF, AFE are similar to each other, and to 
the triangle ABC: and the triangle DEF is the triangle of least 
perimeter, which can be inscribed in the triangle ABC. 

We shall first prove that two adjacent sides of the triangle DEF 
make equal angles with the side of the triangle ABC: that is, that 
BDF = EDC, DFB = EFA, AEF = DEC. 

We have BDF = DFC + DCF, 
as it is the exterior and opposite angle of the triangle DFC. For 
a similar reason EDC = DBE + DEB. 

Now, since the three perpendiculars AD, BE, CF pass through 
the same point G, and in the quadrilateral BFGD the angles at 
D and F are together equal to two right angles, the four points 
B, F, G, D are in the cireumference of a circle: and the angles 
DBG and DFG, or DBE and DFC, being angles in the same seg- 
ment, are equal to one another. In the same way we see that 
DCF = DEB, consequently BDF = CDE, and in a similar manner 
we might prove that BFD = EFA and FEA = DEC. 

We shall next prove that EDC = BAC, DEC = ABC, and 
BFD = BCA. For we have, as before, 

EDC = DBE + DEB. 

Now, since ADB and AEB are right angles, they are angles in 
a semicircle, so that A, B, D, E are points in the circumference of 
a circle, and DBE and DAE being angles in the same segment of 
a circle are equal, and also DAB and DEB for the same reason. 
But BAD + DAE = BAE, so that EDC = BAC; and similarly 
for the others. Comparing the triangles DEC and ABC, we see 
that the angle at C is common, the angles CDE = BAC and the 
angle DEC = ABC, therefore they are equiangular and similar. 
In the same way DBF and FAE are similar to ABC, and the three 
triangles DEC, DBF, and FAE, being similar to ABC are similar 
to one another. 

To prove that DEF is the triangle of least perimeter which can 
be inscribed in ABC, we must avail ourselves of a known property, 
that, if from two points without a line we draw two straight lines to 
a point in the line, their sum will be a minimum when the two lines 
make equal angles with the given line. Now, in this case, having 
proved that FD and ED make equal angles with BC, we see that, 
supposing F and E fixed, the sum of FD and DE is the least of all 
those which can be drawn from F and E to any point in BC. The 
same holds true of the other points F and E, so that on the whole 
the perimeter FED is the least of all those which can be inscribed 
in ABC. 

Various other properties may be demonstrated of this triangle 
of least perimeter. If we call the angle BAD =a, EBC =, 
FCA=y, it is easy to see that FED=2a, EFD=28, FDE=2y, 
and therefore 2(a+$+-y) is equal to two right angles, ora+$+y 
is equal to a right angle. From the construction of the figure it 


appears that a =3 — B, B= 4 -C,y= 5 — A. Also, since 
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to the angles at D, E, F are bisected by the lines intersecting at G, 
ad that point is the centre of the circle inscribed in the triangle 
—_ of least perimeter. If we draw perpendiculars from the ver- 
uF tices of the triangle ABC to the sides of DEF, it will be seen 
me from the similarity of the triangles that they will divide the angles 

into the same parts as the perpendiculars on the sides of ABC, but 

in an inverted position. Now, calling the angles as before a, B, y, 
or AD makes angles a and 6 with AF, AE respectively, and therefore 

the perpendicular on FE makes angles a and /3 with AE, AF re- 
h spectively. The position of AD may be determined by expressing 
it the relation between a and 8, which may be put under the form 
s & =9(a). Similarly, the position of BE will be determined by 
S y = $,(G), and the position of CF by a=4¢,(y). But since these 
lines intersect in one point, one of the three equations must be 
4 deducible from the other two, as the combination of two at their 


point of intersection must coincide with each other. Now the per- 
pendiculars on the sides of the triangle of least perimeter being 
defined by the same equations, of which one is derivable from the 


we take the actual relations, which will be found to be 
a=C—B + p 
B=A—C + ¥ 
y= B— A + a, 
the latter of which is clearly derivable from the other two, so that 
they are not independent, and therefore the lines expressed by them 
must all pass through one point. 
The perimeter of the triangle DEF may be easily found. For, 
from the similarity of the triangles ABC, DEF, we have 


FD: BD = AC:: ABB=5B:e, or FD = “BD, 


but BD = ¢ cos B, so that FD = 3 cos B. 
Similarly, DE = ¢ cos C, and FE = a cos A, so that if p be the 
perimeter of DEF, 
p=acosA + bcosB + c cosC. 
We may likewise obtain an expression for the perimeter involving 
only the sides of the triangle ABC. The area of the triangle DEF 
= 4FD.DE sin 2y = } be cos B cos C sin 2A. 


Since y= 4 — Aand 2y = 7 — 2A. 


Now if r be the radius of the circle inscribed in DEF, its area 
= dpr. 
2 
But 7=GDsiny = GD cosA = CD tana cosA 
cos A cos B cos C 
sin B 





= bcosC cosA tana = b 





other two, must all pass through one point. This may be clearer if 
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5 SEO 

































160 Demonstrations of some Properties of a Triangle. 





Therefore 
cos A cos B eos C 


L Ht itech 
2 op sin B 


= 1 be cos B cos C sinQA 
= be cos A cosB cos C sin A, 
or p = 2esin A sin B. 


Or, substituting for sin A and sin B, their values in terms of the 
sides of ABC, 


1 © mies = ; 
p == po $S(S—a)(S—b)(S—e)} 


8 re . . 
or P — = sS(S—a)(S —b) (S—e)}. 


If R’ be the radius of the circle described round ABC, we have, 
by taking its value in terms of the sides of the triangle, 
abe 
gk" 

Also, if R be the radius of the circle inscribed in ABC, we know 
that generally, 


P= 


abe _ abe 
at+tb+e ~ _? 
if P be the perimeter of ABC. Substituting the value of R’?, from 
this we have 


2RR’ = 


were p2 R? 
tenes abe * 
And therefore if K be the area of ABC, 
8K2 
P = ale’ 


The area of DEF is very easily found in terms of the area of ABC, 
for we found before that it was equal to 


4 be sin 2A cos B cosC = be sin A cos A cos B cos C 
= 2K cos A cos B cos C = k suppose. 
If a’, b', c’ be the sides of DEF, 
a =acosA, b =bcosB, ¢ =c cosC, 

and therefore substituting for cos A cos B cos C, 

ab'c 

‘abe * 

But we have abe = 2PRR’' = 4KR’, 


k= 2K 


abe 

“QR 
It will be seen at once that the areas of the triangles FAE, DBF, 

DCE, are respectively 

K cos? A, K cos? B, K cos?C; 


and therefore = 
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as the sum of these together with DEF make up ABC, we have, 
dividing by K, 
f ab'e 
cos? A + cos? B +- cos? c + g = 
abe 
A singular relation exists between the radii of the circles described 
round ABC and DEF. Let them be R’ and’ respectively. Then, 
a = QR’ sin A, a’ = 2r' sin 2y = 27’ sin QA. 


Therefore a’ =acos A = 47’ sinA cos A. 
Whence t =: Or’. 
We found for r the radius of the circle inscribed in DEL. 
b 
=; cos A cos B cos C, 
sin B 
a b c ’ 
and as ——- = .— = —, = 2h’; 
sin A sin B sin ¢ 
— pee = me 
therefore 7 = 2R’ cos A cos B cos C — R' —. 
K 
, ok 1K 
Also, since r= — = —cosA cos B cosC, 
P Pp 
2K K 
R’ = —, andr = .-. 
P P 


Since the sides of any one of the triangles round DEF are equal 
to the sides of the original triangle multiplied by the cosine of the 
corresponding angle, all similarly situate lines in the triangles, and 
therefore the radii of the circumscribing circles, will have the same 
ratio. 

Now the circle described round F AE passes through G, as A, F, 
E, G, are in the circumference of the same circle: and as AEG, 
AFG are right angles, AG must be the diameter of the circle, 
therefore AG = 2h’ cos A. 

Similarly, BG = 2R’ cos B, CG = 2R' cos C. 

We find also, 

GD =2R' cosB cosC, GE = 2R cos A cosC, GF = 2R’cos A cos B. 

To these geometrical properties we may add the following optical 
one. 

If the interior sides of the triangle be reflecting surfaces, and 
a ray of light be incident on one of the sides in a direction paralle| 
to one of the sides of the inscribed triangle as to FD, it will, after 
two reflexions from each side, return to its original direction so as 
to continue the same course ad infinitum. This property depends 
of course on the equality of the angles of incidence and reflection, 
so that the ray being incident parallel to one of the sides of FDE 
will continue its course successively parallel to the other sides. Let 
fd be the incident ray which being incident parallel to FD will be 
reflected at d parallel to DE, and falling on AC at e is reflected 
\ 
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parallel to EF meeting AB in f’. Proceeding from that point in 
the same manner, it will after reflection again from BC and AC 
meet AB again in some point f”: our object is to show that f” 
coincides with f. For if it does so, it is clear from the course of the 
ray being parallel to the sides of EFD that it will go over the same 
track again, and so on. 
Since the triangles Bdf,; ABC are similar, we have BC.Bd = AB 
Bf. For the like reason we have 
BC .Cd = AC. Ce 
AB. Af’ = AC. Ae 
AB.Bf’ = BC. Bd' 
AC.Ceée’ = BC. Cd 
AC: Ae’ = AB. Af”, 
Adding these together, and observing that we have 
Bd + Cd = BC, Ce + Ae = AC, Af’ + Bf’ = AB 
Ae’ + Ce’ = AC, bd’ 4 Cd’ = BC, 


we have 
BC? + AB? + AC? = AC? + BC? + AB. Bf + AB. Af”, 
whence AB? = AB(Rf + Af”), 
or AB = Bf + Af”. 
Therefore Af= Af” and f" and f” coincide, and the proposition 
is manifest. 
This and several of the other properties of the triangle proved 
here are due to Professor Wallace of Edinburgh. 
V. 


IV.—ON THE SOLUTION OF CERTAIN PARTIAL 
DIFFERENTIAL EQUATIONS.* 


INTEGRATION of the Partial Differential Equation 
1Z\™ dZ\" 
pe ) 

xeE"2 (x) (Fy A= 0, 


a, b, c, m, n being any quantities whatever. 





This equation may be put under the form 


( —eigpetn GLY" [2 2. dz)" 
m Jmtn —_ e a 7Jmtn — —- A = 5 
Xm Zen iy «YE ZMH 0 





* From a Correspondent. 
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m—a 


x™ 


m 





Let x = xa = 
et x= [X-™gX oe mone 
b 


y=f/Y "dy = 


n—b c mtn+e 
= <r m+n 
sd yy z = {ze dZ = m+n 7, mtn. 
n—b : M+n+e 
and the equation becomes 


dz\™ (dz\" 
= — —_ Oise SD 
(5) (5) ‘. Q) 


which is easily integrable. 

If m = a, then x = log, X; if nm = b, then y = log, Y; and if 
m+n -+c= 0, then Z = log, Z. 

But if m-+- 2 = 0, and c¢ does not = 0, then (and in this case 
only) the above substitution for Z will not answer, and equation 
(1) must stand thus, 


aZ\™ (dZ\" 
= ata i) 1D ccceeud 9). 
(=) G) ” (2) 


But by considering x as a function of the independent variables 
Z and y, we may transform this equation so as to be easily in- 
tegrated. 








dx d: 
For dx = < dZ + — dy ; 


dZ dy 
dx 
Ge. = = — dy ap: 
7] dx = dy ; 
dx dZ 
dZ, dx dZ dx dx 
4 nee —- = 1—+-—, — = SS 
and hence = - and > % °7 


and equation (2) becomes by substitution (since m + » = 0) 
da\" 
c — n— 
(5) Zz —A(—1)"=0, 


the integral of which is 
1 


:=s— yZ ™A® + F(Z). 
The same process may be applied when there is any number of 
independent variables. 


dz dz 
Ex. 1. S = a ay*2', 
2 d: dz\2 
This may be written 7? ay = (22 =) : 
3 ‘4 . 
put . =, Qn = 2’, and = = 2’, and we have 


dz dz? 
dy’ ~ dx?” 






oF etree ne: 


a 


f 
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To integrate this we proceed as follows: 


dz' = =f da’ + - dy = pda + qdy' suppose ; 
‘2 = pa + qy' —/(a'dp+y‘dq) = pa' + p’y’ — S(@ + 2py’) dp; 
*, a 4 Qpy' = F(p), 1° 
and 2 = px’ + py’ — F(p), J 


from which p must be eliminated. 


If we restore the original variables, these equations are 


2Vx + gy'p = F(p), \ 
2A = 2V ap + dy*p? — F(p).J 
dz dz 


ix. 2. Zs 


dy adeal y 


this comes under the case where m + » = 0, and c¢ does not = 0. 


Put > = w', and the equation becomes 
: dz dz 
dy = dx’ 
dx’ yal dz + ~ dy ; 
dz — 
= (1 <) dx’ — (= ~) dy; 
dz dy dz 
: a =j]+ e. and = — = = = 
da’ ” dz dy — dy ~ dz’ 


and our equation 








V._ON THE VARIATION OF THE LONGITUDE 
PERIHELION IN THE PLANETARY THEORY. 





OF 

ip; Bie asf 2 25 ae 

Tue following is a short method of finding the expression for 7 
cf 


The usual notation is employed. 


We have the equation 


*= 4 1 +e cos (0 — wf asnabawe (1); 
he 
and since 1—e’= ll 
pa 


e is a function of a and h, therefore 
ea 24 Saye See 
dt do dt da dt dh dt dw dt’ 
dr dr dé 


but — = — —; 
oe dt do dt 


| f 1 lsat dr dr 
therefore, observing tha = — 
eens P 5 dw do’ 


dr dw dr da dr dh 


do dt  dadt dh dt’ 


Now = =_— = —, and = = — = therefore 
dr dw _ 2a* dr d(R)~ dr dR 
“Dhol g! Se: ee a 
d(R dk dr | dR do , 
and = | + 0 i. therefore 
dr dw _ ___2a° dr dr dR 2a* dr dé dr\ dR 
do dt . da dt dr (= da dt + a) do 


le 


p ® da dd dr p 7 da” dh 
dr dd i do _h 


2a° h dr dr dR 2a° h dr dr\ dR : 
—_ — qo 62) 


— —, and —=-. 
dé dt dt 7 
Now, differentiating (1) relatively to a, we have 

1 dr ph de 


tti vend 
utting —= 
I 8 dt 


_ — = — — cos (9 — w).........(3); 
r da hk da ( " (3) 
' Z h? 
and since 1 — e? = : 
pa 
2 . o 
2 ‘ il , and — 2 * = 


= dh 


pa 
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Theory. 
de h de 
whence ae a Qa dh’ *4 (3) becomes 
2a® h dr de 
— - Pm =a == al tae (0 —_ W) oo sererevee(4)e 
Also, r= a(l1 —e cos u), 
“= nt + e—w 4+ esin u; 
d(r) de du 
therefore a ae aa cos u -+ ae sin y — , 
du de du 
and a= Gnu + ecosus, 
wae du __ de sin u 
ho ee —ecosu’ 
d(r) fy sulle de (cos R esin'w \ 
dh dh mar are cos cos u/ 


de cosu—e 


=—2 — —_ 


dh | ~~ # een e 


= —aF cos (0 — w); 
hence, by (4), 
2at h dr d(r) 
pp ® da dh’ 
Substituting in (2) 
dr dw _ dr qk d(r) , dR (oO - a) 
dh dh 





do dt dé dr - dé 
d(r) dr d (9) 
et =5 en a = 


therefore a= “dh dd dh 
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VI.—ON THE OSCILLATIONS OF A FLOATING BODY. 


In all the treatises on Hydrostatics in use in the University, the 
vertical and angular oscillations of floating bodies are treated of in 
separate propositions, on the supposition that either of the two kinds 
of oscillations can take place without the other. It is shewn in 
Webster's Theory of Fluids, Art. 69, that unless the vertical line 
through the centres of gravity of the body and of the fluid displaced 
by it in the position of equilibrium passes through the centre 
of gravity of the plane of floatation, that is, unless the body be 
symmetrical with respect to this line, there cannot be angular 
oscillations unattended with vertical, and it might as readily be 
shewn that there cannot be vertical oscillations unattended with 
angular. It is not difficult, as will be seen, to determine the simul- 
taneous oscillations of both kinds, provided that the angular oscil- 
lations are confined to one plane, namely, that passing through the 
three centres of gravity above mentioned, which may evidently be 
the case when the body is symmetrical with respect to this plane. 

In Fig. 2, AB is the projection of the plane of floatation in equi- 
librium, a6 the plane of floatation at any time, C the centre of 
gravity of the plane AB, a’ b’ parallel to ab through C, G the centre 
of gravity of the body. Let 2 be the vertical distance of the centre 
of gravity of the body below its place of rest at the time ¢, 0 the 
angular deviation from the position of equilibrium. Then the equa- 
tions of motion are 


dx M’ —M 
= — ———- 9, 

dt? uM 7 

ad M'y 

—_—_ — —— q> 

dt? M?? * 


M being the mass of the body, M’ that of the fluid displaced, y the 
horizontal distance of the centre of gravity of the fluid displaced 
from that of the body, MA? the moment of inertia of the body about 
its centre of gravity. 

We must express approximately M’ and M’y in terms of x and 6. 
Let Gh be perpendicular to AB, CA to abd, and let K be the area of 
the section AB, a the distance Ch, a the angle CGA, p the density of 
the fluid. 

Then M’=p x vol. aDb =p x (vol. a’Db’ + vol. aa' bb), 
and vol. a’Db' = vol. ADB, (Webster’s Theory of Fluids, Art. 68) 

therefore p x vol. a Db'= M. 
The vol. aa’ b'b=K x Ch nearly: 


now 2x is the sum of the distances through which the centre of 


gravity is depressed both by the angular and vertical motions. In 
the figure the centre of gravity is raised by the angular motion, 
therefore 

a= Ch — CG cos a — cos (a + 0)3 
= Ch — 0. CG sina, nearly : 
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therefore Ck = x + a0, 
M’=M + pK (x + at), 


dx _—s pkg ; 
aa ee (2 + 29) ....0000 (1). 


Again M'y = moment of aDb= moment a’Dh' + moment a’a’ bd. 
As is proved in all treatises on Hydrostatics in determining the 
equilibrium of a floating body, 

moment a’Db' = (pl — Mb) a, 
I being the moment of inertia of the plane AB about an axis 
through C perpendicular to the plane ADB, and 6 the distance of 
the centre of gravity of the fluid displaced below that of the body. 

The moment of aa b'b = pK x Ch x CG sin (a + 0) nearly 

= pK (a + a?) a nearly, 
therefore M'y = pKax + $p (I + Ka?) — Mb} 0, 
dé yKagq So (1 + Ka?) — MB? 
nee ee D2 ip (I + Ka?) — MB} g 9 (2) 
We Mae Miz baveeessiGeen) 
(1) and (2) are simultaneous equations which may be solved in the 
usual way, and thus a2 and @ expressed in terms of ¢#. It will be 
i Ke I Ka?)—N 
convenient to assume f 7 9 = m?, and p( + a ) Mé =c, then 
a oKa 
they become 


dx 
778 + m? («+ ad) =0, 
¢ 


yn 9 
au mea 


= + ja (@+e0)=0. 


al- 


VIL—NEW ANALYTICAL METHODS OF SOLVING 
CERTAIN GEOMETRICAL PROBLEMS. 


Tue problems which will form the subject of this Article, are some 
in which it is required to show, that the points of intersection of 
certain lines, straight or curved, described according to given con- 
ditions, lie in some other assigned line. The methods of solution 
will be best explained in applying them to particular propositions. 

The first that we shall prove, is one from the Senate-House 
Problems of the present year: namely, 

If on any three chords drawn through the same point in the cir- 
cumference of a circle, as diameters three circles will be described, 
the points of intersection of these circles, two and two, lie in the 
same straight line. 
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Let a be the radius of the principal circle, a, 6, y the angles 
which the chords make with the diameter which meets them; let 
the point through which the chords are drawn be taken for the 
origin, and the diameter of the principal circle for the axis of x. 

The length of the chord whose inclination is a, is 2a cos a, and 
the coordinates of its middle point are a (cos a)*, a cosa sin a: 
hence the equation to the circle described upon it, is 


fa —a (cos a)*2* 4. $y—a cosa sin af? = a® (cos a)’, 


or 2*4y? —2ax (cos a)? —2ay cosa sina = 0...... (1). 
Let tan a, tan (, tan y =), p, v respectively: then (1) becomes 
p? (1 +A?) — 2ax— Qayr = 0.......44 (2), 


putting 2?+4-y?=p* for abbreviation. Similarly the equation to the 
second circle is 


p" (1 +p?)— ax — Daye = 0... .00000(8). 
Subtracting (3) from (2), and dividing by \ — p, we have 
p? (A+ p) —2Zay = 0... ..000 (4). 


Subtracting (1) multiplied by », from (2) multiplied by A, and 
dividing by A — p, 

p? (1—Ap)— Pax = 0...... ..(5). 

Since these last two equations are obtained by combining (2) 
and (3), they will be satistied by the coordinates of the points of 
intersection of the circles represented by those equations. The 
next object is to find a relation between A+p, Ap, and symmetrical 
functions of all the three quantities \, pw, v. For this purpose, let 

B—p?+gi—r=0 

be the equation whose roots are A, p, v: then 
\3— pr?+grA—r = 0, 
we — pee + qu—r = 0. 

Performing on these equations the same processes as before on 

(2) and (3), we obtain 
NMrAwt we —PA+H)+9=0, 
or (A+ p)?—d\p—PQ +4) +.9= 0, 
and Ap (A+ 4)—prAp+r = 0. 

Substituting in these two last equations the values of \ + and 

Ap derived from (4) and (5), we have 


4a? Qax 2ay a 
Pernt wah igen ate 


2ax 
and (“ — ?) (1. -—F)+r=0: 


ms 








+3 @—m)+9—-1=0 vee Wer We (6) 


mugs fp (PRY) — P+ P80... sees (7). 
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These two equations, on account of the manner in which they 
were obtained, must both be satisfied by the coordinates of the 
points of intersection of the circles represented by (2) and (3). 
But, since they involve , », vy symmetrically, they must equally be 
satisfied by the coordinates of the points of intersection of the two 
other pairs of circles. And since their number is the same as that 
of the unknown quantities, x and y, we may, by combining them in 
various ways, obtain every other equation that is satisfied by those 
coordinates, and among them, the required relation. This is done 
by multiplying (6) by 2 and (7) by y and adding; whence, observ- 
ing that 2* + y* = p*, we have 

(q—1)x—(p—r)y+ 2a=0, 
the equation to a straight line. Substituting for p, g, r their values 
in terms of the roots, replacing A, p, v by tan a, tan A, tan y, and 
reducing, we find 

x cos (a tBh+y) + ysin(a+/+y) = 2a cos a cos f, cos y. 

The same method applies to the following problem, of which 
several solutions were given in the Philosophical Magazine in the 
last two years. 

If three tangents be drawn to a parabola, the circle described 
through their points of intersection will also pass through the 
focus. 

Let \, p, v be the tangents of the angles which the tangents 
make with the axis, then their equations are 


" m 
y — Ar = 
y a X 
m 
Yy = px 4 
fA 
m 
y= +-. 
y 
Subtracting the 2nd from the Ist, and dividing by \ — p, 
m 
0=r—--. 
Ape 


Multiplying the 1st and 2nd by A and yp respectively, subtracting, 
and dividing by A—p, 
y = (A+ p)® 
Now if & — pi?+qi—r=0 
be the equation whose roots are A, p, », we have, as in the former 
problem, 
A+ pP—mw—pAtp+9=0, 
and (A+ p—p)\w+r=0. 
Substituting the values of \ + » and Ap in terms of x and y, 
y? — mx — pry + ga? = 0......... (1). 
my — mpx 4 ra? — 
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We must combine these equations so as to obtain the equation 
toa circle. For this purpose, first multiply (2) by y and (1) by 
m, and subtract; then, after dividing by x, 

rey — mgr + m? = 0.........(3). 
Next, multiply (1) by r and (3) by p, and add, in order to 
eliminate the term ay, then 
ry? — rqu? — m(r + pq) « + mp =0.. ......(4). 
Finally, multiply (2) by 1 —gq, and add it to (4), therefore 
r(y? + x) 4+ m(l—q)y—m(ptr)x+ mp =9, 
which is the equation to a circle. Itis easy to see that it is satisfied 
by the values y = 0, x =m, therefore it passes through the focus. 

If a, 2, y be the angles which the tangents make with the axis, 
and for p, g, 7 be substituted their values in terms of tan a, tan {, 
tan y, after putting the equation under the form 

r(y?+22—m?*) +m(1—q)y—m(p +r) (a—m) = 0, 
it becomes 
sina sin f sin y (y?-+2?—m?) + mcos(a+B+y)y 
—™m sin (a+f+y)(*a—m) = 0. 

The succeeding problems are of a somewhat different nature, 
and therefore the same method is not applicable to them, but that 
which we shall employ, is, we believe, in some degree new. The 
first is from the Senate-House papers of this year. 

HKL, PQR are two triangles, prove that if the straight lines 
HP, KQ, LR meet in one point, the intersections of KL, QR; 
LH, RP; HK, PQ lie in a straight line. (Fig. 3.) 

Let O be the point in which HP, KQ, LR meet, and take it for 
the origin of coordinates ; let ab, a'b', a'b", be rectangular coordi- 
nates of H, K, L, and AB, A'B’, A’B" of P, Q, R. 

Then, since O, H, P are in a straight line, 

A B 


tee 
Let this ratio be a, and similarly let 
A’ B ' A’ B 
7=y7=” and al 
The equation to HK is 
(b—U') x— (a—a’)y = ab— al......... (1), 
and that to PQ, 
(B—B’) a — (A—A’)y = A‘B — AB, 
or (ab—ab') x — (aa— aa)y =ad (a'b—ab').........(2). 
Let nx—my = 1.........(3) 
be the equation to a line passing through the intersection of (1) 
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and (2); then, eliminating 2 and y by cross multiplication, we 
obtain the following equation between m and 2, 
$aa'(b—b') —(ab—a'b' 3 m — 3aa'(a—a')—(aa—a'a')in 
4+a—a =0...... ..(4). 
The same line will also pass through the intersections of the other 
pairs of lines, if the same values of m and 2 satisfy these two other 
equations 


$a'q”(U' - b’) _— (a’b! — a’b') m— $a'a" (al —a')—(a'a' —a’a" ys n 


ta’ —a’ = 0......00 (5), 
$a" a(b" -_ b) = (a’b" —ab)} m— $a"a(a" me a)— (a’a" rae aa) n 
40° —a = 0....0500 (6); 


but this is the case, for whether we add (4), (5), (6) as they stand, 
or after multiplying them by a, a’, a” respectively, in both ways we 
get the same relation, whence it follows that any one of the equa- 
tions is a consequence of the other two. To find the value of m, 
multiply by a’a’, aa, a‘a’, and add, then m disappears, and we find 
m=([(a"—a’) aa4(a—a") aa’ +(a' — a) a’a" | 
divided by 
aa a” Sa (U'—b") 4.a'(b" —b) +a" (b— by —(a'b'—a"b") aa 
—(a"b" —ab) aa’ —(ab—a'b’) a’a"}. 
In like manner, 
n= [(a” _ a’) ab b (a _ a") abl! + (a’ oe a) ab") 9 
divided by 
"aala® $a (U = b’) 4 a’ (b" — b) } a’ (6—b'j3 a (a'b' _ a’b") aa 
—_ (a’b’ aaa ab) aa’ —_ (ab — ab’) a’ a 

If the tangents drawn to every two of three unequal circles be 
produced till they meet, the points of intersection will be in a 
straight line. 

9 9 

Let (a—a,)?+(y—B,)” on: " 
be the equation to one of the circles. The equation to a tangent 
line making an angle 0 with the axis of a, is 

(#—a,) sin 0—(y — B,) COS O = yyereceeeee (1). 
If the same line touch another circle whose equation is 
(r7— a,)? +- (y -— [3g)? = Yo" 

we must have 


a, sin 0 —/3, cos 0-4 y, = a, sin 0 — B, COS B+ yoy 


which will give two values of 0. 

To find the point of intersection of the two tangents, we observe 
that if 2 and y in (1) belong to that point, that equation must be 
satisfied by both the values of @ given by (2). Therefore (1) and 
(2) must be identical, whene: 








» we 


° (4). 
ther 
ther 
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wa, __a—a, y—B, _ B,—Bb, 


Se SS > —_— > 
N Yi Ye N Yi— Ye 
a A 
and ¢ = 241%) = Byn—Prir, ‘ 
Viv ys ae 
Let mx 4+ ny = 1 


be the equation to a line passing through this point, then 


(a,y, — 4,7.) M+ (Boy, — By yo) ® = Yy— Yq eee eee eee(B)s 
and the same line will pass through the intersections of the other 
pairs of tangents, if values of m and 2 exist which satisfy also the 
two equations 

(a37q —Ayy,)m +(Byv.—B.%)” = Ya Vyresceveee (4), 

a] 

(a,7,—a57,)™ + (Byy3—Byy,)" = V4 — Yyreeeeeee(5)5 
but this is true, as is seen by multiplying (3), (4), (5) by V9 Yv Ye 
and adding, when the whole will vanish. The value of m is found 
by multiplying by 6,, 8, , and adding; and that of x by multi- 
plying by ay, a,, a, Thus we find the equation to be 


$3, (v2 eee 3) + 3,(y, _ v) +P3(n —y)it 
= 34, (7-73) +4207, — 1) +43(M—Va)hY 
= (a3, —a,,) y, + (498, — 4,33) ¥. + (4,8,—4,8,). 


S. 8. G. 


VIIIL—ON THE INTEGRATION OF SIMULTANEOUS 
DIFFERENTIAL EQUATIONS. 


In the present article we shall apply to Simultaneous Linear Dif- 
ferential Equations the principles, the applications of which we have 
developed in the preceding Numbers of this Journal. The usual 
method for solving these equations was first given by D’Alembert, 
and has received but little improvement since his time, although it 
is so long and tedious, that some change was highly desirable. The 
process which we shall give here is at once simple and direct, and 
shows the advantage of recurring frequently to the principles on 
which our calculations are founded. ‘The theory of the method is 
sufficiently simple. Since we have shown that the symbols of dif- 
ferentiation are subject to the same laws of combination as those of 
number, they may be always treated in the same manner if the 
coefficients be all constants, which is the only case we shall con- 
sider. We have therefore only to separate the symbol of differen- 
tiation from its subject, and then proceed to eliminate one of the 
variables between the given equations, exactly as if the symbol of 
differentiation were an ordinary coefficient. Thus the difficulty of 
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elimination becomes reduced to that between ordinary and algebra- 
ical equations ; and all the facilities for effecting it which have been 
invented for these in particular cases, may be used for differential 
equations. When the equations involve various orders of differen- 
tials of a variable as well as multiples of it, the symbols of differen- 
tiation and of number must be grouped together as one coefficient 
for the variable. Thus, if we had 


d*x dx d*y 1 dy T 

| — —+ce— +m ms 
al ee Mate be bel 
d*y ,dy d*x dx 
_w a }' cd iia 1S , 's _ v : 
dt? ds dt dt? + dt 7a ee T 


they are to be written 
io d ) ‘ pe. d 7 
omen —_ _— c¢-— ae = ’ 
dt ie dt tina + ( dt dt + )y 
d? d ‘ d* ‘ d ' rv 
_— a b - — == > 
and then treated as if each symbol of operation were a factor. 


We shall proceed to give some examples ; and first let us take 
dx + 0 
— ay = 
dt Y : 


dy 


>= 0. 
dt +00 


d —e 
We have to separate = from the variable, and eliminate one 
C 


‘ ae : 
of the variables y or 2, as would be done if a Were an ordinary 
constant. This will be done if we multiply the first equation by 
d . 
— and the second by @, and subtract, when we obtain 
dt 
d* 
— ab) x = 0. 
dt? 

(It is to be observed, that the word “multiply” is used, not 
because the operation is really multiplication, since that is a 
numerical operation, but because it bears a close analogy to mul- 
tiplication, and is represented symbolically in the same manner.) 

Having now eliminated y, we may integrate the equation in a at 
once. The result is 

} 4 
> — ¢ gab)? t —(ab)3t 
r=CyE + Cg ° 


And from the first equation we get 


1 1 
b\} ! b\3 2 
y=c,\— Ps (ab)2t __ ¢, st gia) t. 
‘ 2 a a 
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ebra- As another example, take the equations 
been ds 
ential a 4+ az + by = 0, 
eren- 
eren- dy : ae 
cient an. > oy = 0; 
which may be put under the form 
d 
a\—+a)zx by = 0 
(5 + ) + by , 
d 
a,x — b,\y = 0. 
ye + G + : y 
, ‘ d 
Multiply the first by 7 + 6,, and the second by 4, and subtract ; 
C 
d d } 
diet  ( +b) — able =o, 
1€ {e +a dt + 5, a" | 
The coefficient of a is obviously of the second order. Let it be 
made up of two factors, so that 
d 
a h k) « sz Q 
(5 + ) (5 + , 
where h +k =a+b, hk = ab, — ab; 
so that they are the roots of the equation 
2? — (a + b,)z + ab, — ab = 0. 
Integrating in the usual manner, 
ne w= Ge ht 4 cg et; 
and from the first equation we deduce 
uy h—a k — 
aah. en ee 
by eT, + ae 
It will make no difference in the method of the solution, if 
there be a function of ¢ on the other side of the equation. As, 
for instance, we have 
dx ; . : 
— % Sy = 
ot dt + v + 17] > 
I Y on 4. 5 
= — H i OY = 
| dt + + oY &, 
at or (5. 4) 2 + Sy = ¢ 
d + Yy 9 
(; ao s)y + y+ uw = ef. 
Eliminating y, 
fd /d +5) 
— —6sa=—1 + 5t — 3¢ 
(a + *) (G+ } ‘ins 
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d d ) 
—~+2)(—+4+7 = 5t — Se. 
or (5+ ) (47 x ] + 5t Se 
: ; d is 
Integrating with regard to = + 7, 


« 2) _ 5, e. - 
at oma Tet 2 ae 


$1 5 1 
and z= — — — t— —-¢! ce + oe-%, 
i96 * 14 rere}! 
the value of y will be found to be 
9 ] 5 8 
=— - - ene KOT ee 
¥~o3 9° tat ~ 35 _ 
Take also the equations 
dx a . 
— 5a y = ¢, 
dt + ey , 
dy 
—- Sy — r= 2, 
7, + SY E 


Eliminating y, we obtain 


] 2 
(++ +) 24° — *; 


d —2 
whence x = (= + t) (4c! — &*) 4 (ct + ¢,) 


dt 
4 1 , 
re = 35  — a* + (ef + ¢,) , 


] 7 
ea oe . —4t, 
and y = og 4. 36° + (et + ¢,)¢ 
If there be three simultaneous equations, the same method is to 
be pursued, but the calculations must be necessarily long, as in 
ordinary elimination. We may, however, avail ourselves of the 
method of cross multiplication. If, for instance, we have 


‘ + by + cz = 0, 
2 +dax+cz= 0, 
dz 


77 +axr+by = 0, 


d 
or au a by ~ eee ere (1), 


d 
ax + a! BOE Be] oo5cee ses cvecee (8), 


ID cas kinias uagnss Dy 


, , d 
ot ay ae 





Wi 


al 


Vv 


> —- 


a se 


J 
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we may eliminate y and z by multiplying 
d* d 1 
(1) by a — bc, (2) by be —b ih’ (3) by be’ — c= 


and the result will be 


{2 (S - v) a 5) (0 — oY a = 
tai (ga ~ 86) +a (Be - 85) +o (eG) ie =o, 


3 

or fd — (ab + a’e + be) + abe + a'be | x= 0, 
: ) dé dt J 

which can be integrated by the usual method. The equations for 

determining y and z will, of course, be symmetrical with that for a. 

There will be three arbitrary constants arising from the integra- 
tion of this equation; and it would appear, that as there are other 
two similar equations, from which six more arbitrary constants would 
arise, there would be nine on the whole. But there are really only 
three independent arbitrary constants, as we are able to deduce 
the other two variables y and z from the value of the first, without 
integration, and consequently the arbitrary constants in their ex- 
pressions must be derivable from those in the first integral. 

The same method of integration may of course be applied to any 
number of simultaneous differential equations ; but the difficulty of 
elimination rises rapidly with the number of variables, as in the 
case of ordinary numerical equations. We shall take as an exam- 
ple of four simultaneous equations, those given by Mr. Airy for 
determining the secular variations of the eccentricity and longitude 
of the perihelion, Plan. Theor., p. 123. 

They are of the form 


on sae : 1) 
= —u SGA das des saesc ) 
dt t 2 ( 
0 : j (2) 
SS GU — VU — gl coc ccccecccvees (2), 
, dt . 
d 
O=—u Biel OD Gacccaicsivnce 4h 
dt + | 2 ( ) 
d 
0=bu' — OC =) CBW sas sisievcesencs C4 
‘ dt . (4) 
To reduce them to three simultaneous equations, eliminate v 
d , 
by multiplying (1) by = and (2) by a, and adding. Then 
L : vi A f 


a 2 ‘ d ‘ - 
ap + M1) & — Ua — a0 = Oi tscsheien COP 
Again, eliminate the same variable between (1) and (3), by 
multiplying (1) by 4,, and (3) by a, and adding, so that 
d d , F ' 
b, at +4, a” + (a,b, — ab,) PD cscsensevene CO) 


AA 
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and the equation (4) is 
d 
bw — byl + — vo HO oo. cceceeceeeee (7) 
iss dt 7) 


Eliminating wv and v' by cross multiplication between (5), (6), 
and (7), we obtain an equation which reduces itself to 


dt* 
which we shall have no difficulty in integrating. 


‘ha 12 5 
‘5 + (a, + b,? + 2a,),) a + ab, — a} 0 ot O cevuns (8), 


From its form we see that the operating factor is decomposable 
into two binomial quadratic factors, so that we may put it under 


the form 
d? 2 ad 2 
(5a + ") (at Su = 0, 


—k,?, — k,* being the roots of the equation 
22 + (a? + 5? + 2a,b,) z + a,b, — a,b, = 0, oo. see (9). 
Integrating with respect to the first factor, as in page 28 of our 
first Number, 


d? : 
( 1. h?) u =e, coskw + e, sin ha, 


dt? 
= ¢, cos (hw + a); 
by changing the form of the constants ; whence 


ad? 2) ! ; 
“u = (ss + k, j $c, cos (kw + a,)} + ec, cos (kt + ay). 


The operating factor will only affect the constant of its subject, 
and as that is arbitrary, we may write 
u =e, cos (kw + a) + ¢, cos (kgt + ag). 
To find the value of v, we may observe, that if we had elimi- 
nated wv’ instead of v at first, we should have, instead of (6) and 
(7), the similar equations 


: d 
My ou + 4, a + (a,b, — ay) 0 = 0, 
‘ d 
au — au+ a o= 0. 


Eliminating w between these by multiplying the second by —, 
) iP 
and subtracting, 


‘ha 
a — (a), — a,)} v = (a, + 4,) 


du 


a’ 


or putting for w its value, and effecting the differentiation indicated, 
{ d 
| dé 


— (ab, - a,p,) 


ay 


J 
= —(a, + 5,)§e,h, sin (ke 4 a.) + ¢h, sin (hw 4 a); 











(6), 


(8), 


ible 


der 
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whence, putting a,b, — a,b, = A, 
chi 
=(a,+5,) x 5) sin (h, wa)+ pa sin (z, #ta)}. 
By substituting the actual expressions for 4, and k,, this may 
be simplified. For solving equation (9), we get 
h2=h$a,2 4624 2a,b,—(a, 4-b,) [(a,—5,)? + 4a,),] 2: 
then kP+A =4 (a, +,) fa, +5,— [(a,—4,)? + 4a,b oJ] %. 
But we have also 
k, =} fa, + 6,— [(a, - b,)? + 4a,b,]3 i, 
and slaalianly for his so that 
v =e, sin (hy + a,) + ¢, sin (kx + a,). 
Knowing uw and v, we easily determine w and wv’. For, by equa- 
tion (2), 


/ 
au = a Nad — =v 


dt 
=(a,—h,) ¢, cos (kx + a,) +(a, —hk,) € COs (he + ay) ; 
similarly v’ may be found. 
In this case we have had an example of the integration of simul- 
taneous differential equations, of an order higher than the first. 
We shall take another in 


d*y 

— —az — by=ec 
dt* Y % 
d*z hi 7 
— —az—by=e 
dt® isi 


being two of the equations for determining the circumstances of 
the movement of a floating body in a position nearly of equilibrium, 
the other two being 


do 
dé + no —m=0, 


d’s d’ dp sia 
dt® IW dt? : 


Putting the equations under the form 


ad ‘ hi : 
—-—a\)z—6 c 
dt’ ) Y ’ 

d? 


we can eliminate z by multiplying the first by —- a, and the 
de 


second by a, and adding, when we get 


(¢ : (5 ») alas - 
\de «) de = oe Ss ay = ae — ca, 
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If k,?, k,? be the roots of 
2 — (a 4 b)z + ab — ab’ = 0, 


the equation becomes 


a \ (@ Cy 
(Ge) (Gem at) y= ae = on 


the integral of which is 


‘ ‘ 
ac — ca 
i> rn okt > —k,t » » ewkot 
— hh? + ev + Qe + ey + ce ™, 
9 
ac’ — co 


= Ty » owt » gkot » owk 
= Coe “YE 4 Clet2 + C,€ “54. 
ab — ab! iil t % 4 a s 
The value of z may be found from that of y; and if we know the 
initial circumstances, we can determine the arbitrary constants. 


The same method may be extended to simultaneous partial 
differential equations, according to the principles developed in Ar- 
ticle IX. of our third Number. Take, for instance, the equations 


dz du dz 


—+¢—-+a—-4+ bz =0, 
dy 


dx dx 
lz if if 
c ms t - +a ~ + bu = 0, 
dx dx dy 
which may be written 


d d d 
(+a: +b)z+eLu=0, 
dx dy dx 


7. oe d d 4 ) ie 
c z* + (= +a dy 21 oe = 0. 
, : ‘ d d 
To eliminate uw, multiply the first by € +a—4+ ») » and 
: ” \de dy 
_d he 
the second by ¢ dn’ and subtract. Then 
x 
d 1 2) 
{ 2 (a -{- 0) (. — of 0) 
ae dy d 4 dy — 
dx® l—cc dx l—ce sirens 


If we put 1 — Vecd= m,, 1+ Vec= m,, this may be re- 
solved into factors, 


= il n («i 5 + Dies ' m (as * i) } == 05 


the integral of whic hi is 


d d 
i a . dy + ») (y) 4 g-Mat # dy ») yp (y)s 
or 2 = & "7 o(y — am x) + emer Y(y — am,x). 


. Airy, at page 279 note of the Undulatory Theory, gives as 








nd 
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the equations for determining the small disturbances of an elastic 
medium in three dimensions, 


dus, dd (du dv dw 
de" dx G dy =) 
d*v_—s «ss hs (dus dv dw 
ade ~ . dy (= dy 4 


d°*w _d (du dv dw 

— =a —|— ca as 
; dt* dz \dx © dy ~ dz 

We might eliminate v and w, by cross multiplication, between 

these equations, and so obtain an equation in w which might be in- 

tegrated, but it will be more convenient to proceed as follows. Let 


—— ( + dv + =) 
ou dx “ dy * dz]’ 
d" du d— dv d™ dw 
then! Ss ee ee gs. 
dx" d@ dy" d@ dz" dé 


ae & a de 

ax? dy?’ i’ and add; then 

d* (du dv dw ad a? da? 
de es a 4 - (a bs ay + *) 


ad ¢ . ¢ \ 
ey ee , 
or ie a (a + -— dy? + dz mir 7 


the integral of — is 


Multiply by 


- 
1 (Sta ay* dz a p (a, Y z) + € coat (E+ apt dz x) W(a, 4 Ys 2 
knowing 7, we can determine u, v, and w, since 
ad dr d~* dr d-? dr 


adi dt=2 dx’ ps dt~ dy’ _o dt-? dz 


This solution is due to Mr. Greatheed. 


We have now given a sufficient number of examples to enable 
the student to understand thoroughly the method, and we think 
that they show clearly the advantages of a process, which, to some 
persons, might appear to carry out to a startling extent the princi- 
ples on which it is founded. 


D. F.-G. 














IX.—ANALYTICAL SOLUTIONS OF SOME PROBLEMS 
IN PLANE ASTRONOMY. 


WueEN it has been required to determine the ratio of the small 
variations of ares or angles on the surface of a sphere, mathe- 
maticians have generally employed a geometrical method, probably 
because it ‘appeared to them that the application of the Differential 
Calculus would lead to laborious processes. But it will be shewn 
in this Article, that Analytical solutions of such problems are not 
necessarily tedious, and the method by which they are abridged 
will be seen to be very similar in all. 

1. When the latitude and hour-angle are determined from two 
altitudes of the Sun, and the time between, to find the errors caused 
by given small errors in the observed altitudes. 

Let / be the latitude, z one of the observed zenith distances, 
h the hour-angle from noon, 6 the Sun’s declination, a his south 
azimuth. Then 

cos z = sin sin 6 + cos / cos 6 cos h. 
Differentiate, using A to denote small finite variations, which will 
be nearly in the ratio of the differentials; and considering ¢ con- 
stant, as the Sun’s declination is supposed to be known accurately : 
therefore 
—sin zAz=(cos/ sind—sin/cosé cosh) Al coslcosdsinh Ah. 
Now cos Z/sin 6 — sin 7 cos 6 cos h =— cos @ sin 2, 

and cos 6 sin A = sin z sin a; 
substituting and dividing by sin z, 

Az = cos aAl + cos / sin adh. 


If 2, a be the zenith distance and azimuth at the second 
observation, 
Az’ = cos a’Al 4 cos /sin a’Ah; 
Az sin a’ — Az sina 
therefore Al = . ’ 
sin (a’ — a) 
Az’ cos a — Az cosa 
and Ah = —> 
cos d sin (a — a) 
2. To find the errors caused by parallax in the hour-angle and 
declination of a heavenly body. 
The notation remaining the same, we have the equation 
cot h sina = cos / cot z + sin 1 cos a, 
in which / and z are the only variables. Differentiating, 


sin a cos ¢ 


‘= 7 , mse 
(sin 2) 
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But if P be the horizontal parallax, 
Az = P sing, 
4 cos 6 sin h 
and sina = ———__; 
sin z 


s Usin h 
therefore Ah = —, in h 





cos 6 
. sin 06 — sin Z cos z 
Again, — cosa = — apie. 
cos Z sin z 
or cos a cos / = sin / cot z — sin 6 cosec z. 


Differentiating, 
O= $sin 6 cosee z cot z—sin 1 (cosec z)*? Az—ecos 0 cosee 2 Ac, 
ji 
0 = P (sin 6 cos z — sin 7) — cos 6 Ad, 
but cos z = sinJsin d + cos / cos 6 cosh; 
therefore cos 6 AC=P §sin Z (sin 8)? 4 cos J sin 6 cos 6 cosh—sin 2, 
and Ad = P (cos /sin 6 cos h) — sin I cos 0). 
3. To find the precession of a star in right ascension and 
. . J © 
declination. 
Let a, 3, J, X, be the right ascension, declination, longitude, 
and latitude of the star, and w the obliquity of the ecliptic. 


Then sin 6 = cos w sin + sin w cos X sin J; 
therefore cos 6 Ad = sin w cos A cos J Al. 
3ut cos \ cos / = cos 6 cos a; 
therefore Ad = sin w cos a Al, 


= 50",18.é sin w cos a, 
¢ being the interval expressed in years. 
Again, sin \ = cos w sin 6 — sin w cos 6 sin a; 
therefore 
0 ==(cos w cos 64sin w sin 6 sin @) Ad—sin w cos 6 cos a Aa. 
Substituting the value of Ad, 
Aa = 50",18.t. (cos w + sin w tan 6 sin @). 
4. To find the aberration of a star in longitude and latitude. 
Let S (fig. 4.) be a star, O a point in the ecliptic 90° behind 
the Earth’s plane, SL perpendicular to the ecliptic. 
Let /, X be the latitude and longitude of the star. Then 
YO = 180° + © — 90° = 90° + ©, 
LO= 90° + © —1. 
By Napier’s rules, 
cos O = cot OS tan LO, 
= — cot OS cot (© — J). 
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Taking the logarithmic differentials, 
as A.OS Al 
~ sin OS cos OS © sin (® — /) cos (© ey l)° 
But A.OS =n sin OS; n denoting the ratio of the Earth’s ve- 
locity to that of light, and 


cos OS = cos SL.cos LO = — cos d sin(@ — /). 
L08 _ 7] 
Hence Al = — n = (O = ; 
cos X 
Again, sin \ = sin O.sin OS; 


therefore cot \ AA = cot OS.A.OS, 
= n cos OS, 
=— ncosAsin(® — J); 
therefore AX = — nsinAsin(® — 2). 
5. The aberration in right ascension and declination may be 
found in the following manner. 
Since tan a = cos w tan / — sin w see / tan \, 
(see a)? Aa=cos w (sec 1)? Al—sin w see 7 tan Z tan XA 
+sec 1 (sec \)* AX? 
—n cos w (sec 1)? sec X cos (© —2) 
+” sin w sec / sec \ tan d jtan Z cos (© —/)+4sin (© —l)}, 
=—™n cos w (sec 1)? see X cos (© —l) 
+ sin w (sec 7)? sec \ tan A sin ©, 
cos w COS © +(cos w tan /--sin w sec J tan A) sin © 
: ls paaaasacatel cos J cos A 7 ' 
cos w cos © +tan asin ©- 


COS @ COs é 


-=— 





; COS w COS a COs © +sin a sin © 
therefore Aa = — n ae ween. 
cos Cc 
Also, since sin d=sin w sin Z cos X + cos w sin X, 

cos 6 Ad = sin w (cos / cos \ Al—sin J sin AX) + cos w cos X Ad, 
=—zn sin w cos Z cos (© —/)—sin 1 (sin X)? sin (© —I) 

—n C08 w sin X cos d sin (© —2)}, 
=--n SIN w COS © 

—n (sin w sin / cos d 4-cos w sin d) cos d sin (© —1), 

——n sin w cos © —” sin 6 cos J sin (O —l), 
=—n sin w cos © —n sin 6 cos / cos d (sin © —cos © tan 2), 
—=—n sin w Cos © 
—nsin 6 cosa cos d §sin © — cos © (cos w tan a+4sin w tand sec ay, 
— --n sin w (cos 0)? cos © 


-n sin 6 cos ¢ (cos a sin © —cos w sin a cos ©); 
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therefore Aa=—n sin w cos 6 cos © 
4sin 5 (cos a sin © —cos w sin a cos © )}- 
6. To find the time of year at which the Sun passes from one 
given zenith distance to another in the shortest or longest time. 
Let Z (fig 5.) be the zenith, P the pole of the equator, S, S’ 
the two positions of the Sun. Let z, z be the given zenith 
distances, h, h’ the corresponding hour angles, 6 the Sun's ¢decli- 
nation, and / the latitude of the place. Then 
cos z = sin/sin 6 4 cos cos 6 cosh...... .. «00 (1); 
cos 2 = sin/7 sind 4 cos / cos 6 cos h’..........00008 (2); 
and f’ — A is a maximum or minimum ; therefore 
0=(sin 7 cos 6—cos / sin 6 cos h) dé—cos 1 cos é sin h dh, 
0=(sin J cos §—cos Z sin 6 cos h’) di—cos 1 cos 6 sin h' di’, 
dh’ — dh = 0; 
sin 2 cos 6 — cos J sin 6 cos h’ 
hence ——s/ 
sin h 


sin h 


sin 7 cos 6 — cosZsin dé cosh\ .. 

G3 =z O 2.020850. (8); 
therefore, either di = 0, which gives the times of the solstices, or 
the other factor is equal to zero. Now 


sin 7 cos 6 — cos 7 sin 6 cosh " 
—— = cot PSZ; 


cos Z sin h 
therefore cot PSZ = cot PS'Z, 
therefore cos PSZ = + cos PS’Z, 

sin 7 — cos z sin 0 sin 2 — eos z' sin 3 
or ; = =a ae — 
sin z cos 6 sin z cos 6 
- . Oa ie : 8 ae : 
sin / (sin 2’ sin z) = sin 6 sin (2’ F z)............ (4). 
Taking the upper sign, we have 


’ 


~ al ~ i ' 
2 








Pe idl , ss LD 
2 sin 1 cos = sin > = 2 sin 6 sin —5— cos aes 
»' + » 
N68 
bg ‘ 9 
therefore sin 6 = sin / . sonsenvonnee CODE 
cos —5 


whence the time of year is known. 
If it be required to find the time between the zenith distances, 
we have from (3) 
sin 2 cos 6 (sin ’ — sin h) = cos / sin é sin (h' — h); 
hi +h Km—h 
¢€ 


therefore sin 7 cos é cos —— cos 7 sin 6 cos as 


BB 
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By adding (1) and (2), 
2’+2 Z—z +h Med 


cos == sin / sin 6+cos / cos 6 cos —-— co 


2 2 
Ath 
Eliminating cos — 





cos 





2 








, 


eq 


U , 


= ~ SS —_ 
“ - ~ ~ 





sin 7 cos : : 
2 m! 


h’ — hy? 
= sin 6 {sin 1)? + (cos 2)? cos Fs) \. : 


Substituting the value of sin 6 from (5), 


a! — 2 h! _ h 2 t\ 
(cos = ) = (sin 2)? 4 (cos 1)? cos (—) ; “ 
,_ #—2% ; hi — hie r 
therefore (sin = ) = (cos 1)? (sin >) . 


_ k—As _— Isi 2 —2Z (6) 
sin —5— = + see 1 sin —Z— o.eereeeeeee (6). 


~ 





To find the results of taking the lower sign in (4), we have only 
to change the sign of z, which gives 
/ 





— : 4 : 
sin 6 = sin 7 : ovesinevensene (4 Do 
"2 +2 
cos — 
2 
_hk—h . £+2 
sin —— =-+ sec /sin QT ote (8). 


Since there are two periods in the year at which the Sun has 
the same declination, within the same half year, and at equal dis- 
tances from the equinoxes, equations (5) and (7) give each two 
periods of maximum or minimum duration, and the solstices are 
two others, so that there are six in all. But one or both of 
the equations (5) and (7) may give a value of é greater than 
the obliquity of the ecliptic, which, therefore, the Sun never 
attains, so that there may be only four or two maxima and minima. 
In the problem of finding the time of longest or shortest twilight, 
the values of z and z’ are 90° and 108°; so that equations (5) and 
(7) become respectively 

sin 6 = — sin / tan 9° ...........0000006 (9), 
sin 6 = — sin J cot 9° ................ (10). 


Of these, the former gives a value of 6 less than the obliquity, 
whatever be the latitude. The value for north latitude, 52°, is 
— 7°.10'.10", which occurs on the 2nd or 3rd of March, and the 
11th or 12th of October. But (10) will give a value of é which is 
never attained, unless 7 be not greater than 3°.36'.55". At the 
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equator the values of 3 given by (9) and (10) coincide, being both 
h zero. 

After some considerations, which enable us to distinguish between 
maximum and minimum values, we obtain the following results. 

At the equator, the duration of twilight is a minimum at the 
equinoxes, and a maximum at the solstices. 

Between the equator and latitude 3°.36.55", the duration is a 
minimum at two periods near the equinoxes, in the winter half of 
a the year, a maximum at two others nearer the winter solstice, a 
minimum at the winter, and a maximum at the summer solstice. 

At and above latitude 3°.36'.55", the duration is a minimum at 
two periods near the equinoxes, in the winter half of the year, and 
a maximum at the solstices. 

The seasons are here spoken of, taking into account their being 
opposite on different sides of the equator. 

C. G. 


X.—ON CERTAIN CASES OF CONSECUTIVE 
SURFACES.* 


THERE is a very simple theorem, which does not seem to have 
been much noticed by English writers, but which is of frequent 
use in Analytical Geometry, particularly when all the equations are 
put in a symmetrical form. It is, that the formula 


a _ a’ a’ 


p“7"s* 
always implies the truth of the following, viz. 
_atadza’t+.,..  matmd + mid +... 
~~ mb + mo + mb" +... 


a’ 
fteeiee 8 FYE 





m8 





Ve pa 4 a? 4 
Vie Oy. 
m, m', m",... being any quantities whatever. The proof is so 
obvious, that it is unnecessary to put it down here. This being 
premised, suppose the equation to a surface is 
u=h oO errerrerrerr err ert ery (1), 


where / is an absolute constant, and wu =F (x, y 2, a, b, c), 
a, 6, e being parameters, by the variation of which the surface 








* From a Correspondent. 
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takes its different positions. And first, suppose a, 6, and ¢ to vary, 
subject to the single condition 
Jes B 0+) eee B ..0.250.00000800s (8), 
which is equivalent to supposing the equation (1) to contain two 
arbitrary parameters. Now, suppose the partial differential co- 
efficients of uw and f(a, b,c), taken with regard to a, b,c, are 
P, Q, R, A, B, C, respectively. Then, in order to find the 
intersection of the surface with its consecutives, we shall obviously 
have to combine the equation (1) with the two following, viz. 
Pda + Qdb 4+ Rde = 0, 
Ada + Bdb 4+ Cde = 0. 
If between these two equations we eliminate da, we obtain 
(AQ — BP) db + (AR — PC) de = 0; 
and since there is no relation between db and da, their coefficients 
in this equation must separately be equal to nothing; that is, 
AQ — BP=0, AR — PC = 0, 
whence we immediately derive 
P Q_R 
This formula, combined with (1), determines the coordinates of 
the point in which the surface is intersected by all its consecutives ; 
and if we also combine it with (2), and eliminate a, }, and e, 
we shall obtain the equation to the envelope. This elimination 
will often be easy, if the functions w and f(a, , ¢) are homo- 
geneous with respect to a, 6, and ec. For suppose m and x are 
their respective degrees of homogeneity, and we shall have 
aP + bQ + cR = mh, 
aA + bB + cC = nk. 
Hence, applying the theorem above laid down to the formula (3), 
we get 
r. © (= aP + bQ 4+ cR\ _ mh 
5 a eateat oo (4); 


' 


A BC 
by the help of which a, 6, and e¢ may be eliminated from either of 
the equations (1), (2). 

If the generating surface was a plane, suppose its equation was 
ax + by + cz = f*...... .....00. (5), 
and we should have P=a, Q=y, R =z; and hence, sub- 
stituting in the formula (3) and applying the same theorem, we 
should get 
a a h? _ Var + y® 4 22 (6 
AB C aA+bB4¢cC Vata BY 
These equations determine immediately the coordinates x, y, z 
of the point in which the plane is intersected by all its consecu- 


|& 
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tives, and the length of the radius vector r from the origin to this 
point. Also, if A, », v are the angles which r makes with the 
three axes, we have evidently 
cosA cosh COSY __ 1 o 
ay — . BR = C a V4 By @ ecce eee (7). 
Hence it is easily seen, that if we described a surface whose 


equation was 
I (ar Ho B) Re... c00000 000000 (8), 
a, b, and e would always be equal to the coordinates of some 
point in it; and if a tangent plane was drawn at this point, the 
perpendicular upon it from the origin would make precisely the 
same angles with the axes as those indicated by the formula (7), 
and would therefore coincide with the radius vector r. Moreover, 
if p was the length of the perpendicular, we should have 
aA + bB + eC 
~ VA? 4 Be 4 C2’ 
and therefore it is evident, from the equations (6), that rp = h'*. 
Hence, if in the perpendicular from the origin upon the tangent 
plane to the surface (3), we take a point, such that its distance 
from the origin multiplied by the length of the perpendicular shall 
= h’, the locus of this point will be the envelope of the system of 
planes defined by the two equations 
ax + by + cz =h’*, f(a, b,c) =k. 
As another example, suppose the equation to the surface is 
9 :. © ¢ 
where ¢, x, and y are any functions of x, y, z, not containing 
a, b, ec; and suppose the relation between a, 6, and ¢ is 


(*)' + (F)" + (2) a ee 


here the equations (4) become 


fe CS ae l : 

- = —+—_ = me SS Dt cceccocee (11)3 
a®q"-n p" bm n y*e" n 

and if we take the values of a, b, c from these, and substitute them 

in (9) or (10), we get the equation to the envelope, viz. 


n 


(s)"" + (x)"" ~ (=) ale 
a!" 3m y 


In the particular case when »=™m, the equations (11) are 
reduced to the three ¢=a™, x = pb", ) = y", and the envelope 
consists of three distinct surfaces defined by these equations. And 
if 9, x, W each contain only one of the variables 2, y, z, the 
envelope is reduced to a certain definite number of points. For 
instance, all the ellipsoids defined by the two equations 
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n2 2 a 2 2 2 
ceed, 6 oe 4 wi, 
a b c a* b° c 
will intersect one another in the eight points, whose coordinates 
are ta, +f, + y. 
But if a, b, ¢ were subject to two equations of condition, 
suppose, for instance, 
F(a, & 0+) 2 Bo... 000000000000 (18), 
J, Ca G, ©) wt By, 000.000000000008 (18), 
(which is the same thing as supposing the equation to the surface 
to contain only one arbitrary parameter,) let A,, B,, C,, be the 
partial differential coefficients of f, (a, 6, c), the others remaining 
as before, and we shall have the three equations 
Pda + Qdb + Rde = 0, 
Ada + Bdb + Cde = 0, 
A,da + B,db + C,de = 0. 
To eliminate the differentials, apply the method of indeterminate 
multipliers, and we shall have 
P=rAA+ rAy) 
Q = AB + pB,)} ............... (14). 
R =aAC + uC, 
Let aP4+bQ4+cR=v, aA+bB4cC=S, aA,+6B,4+cC,=S,; 
then, if we add the equations (14), after having respectively multi- 
plied them by a, 6, and c, we get 
v=)AS + pS); 
and if we substitute in the equations (14) the value of \ derived 
from this equation, we obtain 
SP — Av = »(A,S — AS,), 
SQ — Bu = » (BS — BS,), 
SR — Cv = » (CS — CS,), 


ll 


and therefore 
SP— Av _ SQ - Be _ SR — Cv (15 
A'S — AS, BS — BS, C,S — CS, "°°" > 

This formula, combined with (1), comprehends the equations to 
the curve in which the surface is intersected by its consecutive; and 
if we eliminate a, b, ¢ by means of (12) and (1%), we shall obtain 
the equation to the envelope. If the functions u, f(a, b, c), 
JS, (@ 4, ce), are homogeneous, and of the degrees n, m, m, respec- 
tively, then v = nh, S=mk, and S'= mA. 

Suppose, for example, it was required to find the surface gene- 
rated by the consecutive intersections of tangent planes drawn to 
the ellipsoid, whose equation is 


at, ee 16 
wet GH eS Desc eeeees (16), 


9 
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at all the points where it is cut by the plane, whose equation is 
ke + my + nz ah ........0006 (17) 
The equation to a tangent plane at the point (2, y, z), is 
Le 
ae + RB + > fe DE diedi cove OE 
t,, being the coordinates of any point in the tangent plane. 
Here the variable parameters are 2, y, z, subject to the two equa- 
E 22 
tions (16) and (17). And we have P = —4 A= =? A, =4 
&e., v= 1, S= 2, S, = kh; hence the formula (15) becomes 
E—-rv% — mn-y _ f-2 


ay vereeeee (19), 


which is equivalent to the equations to the generating line of the 
envelope. We may put it under another form, by subtracting 1 
from each term, after having multiplied by ; we thus obtain 

ki — la? — kn — mb? — hG—- ne* 

la? —kx ~~ mb? —hy ne? — kz" 


la? — hx ~ mb? —hy ~~ ne? — 


This shews that the generating line always passes through a fixed 
: : la® mb? ne? 
point, whose coordinates are . ae eo 
k 


In order to eliminate 2, y, z from the equations (19), if we mul- 


tiply the numerator and denominator of each term by and 


6 ) 
a2’ §’ 
4 , 
= respectively, and then add the numerators and denominators, we 
c 
obtain an expression which, by the help of equations (17) and (18), 
immediately reduces to the following, viz: 

2 2 72 

S U] - 

—~+-+--1 

a bi b? ad ce 

E+ mn + 2 —h- 

Again, if we multiply by /, m, 2, and add, we obtain 
“E+ mn + ni —hk 


Pa? 4 mb? 4- n%e? — h*’ 





and these two expressions must be equal, by the principle stated at 
the beginning. Hence, 
r2 
“ 


(lé + mn + nz — hy 


T,¢ 
+> wt hie 
ae + be + c2 Par + mb? + nrc? — k2’ 





which is the equation to the conical surface required. 

If it was required to find the locus of the vertex of this cone, 
when the plane which cuts the ellipsoid is moved, but made to pass 
through a fixed point (a, (6, y), suppose a, y, z are the coordinates 
of the vertex, and we have seen above that we shall have 
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“=! = im “= ere. 
la me ne k 
also, a, 3, y must satisfy the equation to the cutting plane, that is, 
la + mB + ny =k; 
multiplying, then, the terms of this last equation by the terms of 
(20), we have 
o+ e+ f=, 
which is the equation to the locus sought, and represents a plane. 
If the point (a, 3, y) is a point in the surface of the ellipsoid, this 
equation coincides with the equation to the tangent plane at that 
point, as it evidently ought: 


a h? 


M. N. N. 


XIL.—MATHEMATICAL NOTE. 


Geometrical Theorem. Let A,A,A,... be a polygon of x sides 
inscribed in a circle, a,, a, a,, &c. the angles which the sides 
A,A,, A,A,, &e. subtend at the centre. Then 


iz: A,A,A, =a — AA,A,, 


a, + ay 


a 
‘o 


> 
a, + we 
9 ’ 


&e. 


. 


Q,-) + Uy 
2 
If be even, adding all these together, we get 
4 A,A,A, + 2 A,A,A, + &e. + 2 Ax A,A, 
n Qn , . 7 
— — «== _— C _— 
=5tFI-pZv= lm 2) 5° 
or the sum of the alternate angles is equal to » —2 right angles, 
a curious extension of Euclid, 111. 22. 
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